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ABSTRACT 

We present results from X-ray observations of the quadruply lensed quasar 
RX J0911. 4+0551 using data obtained with the Advanced CCD Imaging Spec- 
trometer (ACIS) on board the Chandra X-ray Observatory. The 29 ks observation 
detects a total of ~ 404 X-ray photons (0.3 to 7.0 keV) from the four images of 
the lensed quasar. Deconvolution of the aspect corrected data resolves all four 
lensed images, with relative positions in good agreement with optical measure- 
ments. When compared to contemporaneous optical data, one of the lensed 
images (component A3) is dimmer by a factor of ~6 in X-rays with respect to 
the 2 brighter images (components Al and A2). Spectral fitting for the combined 
images shows significant intrinsic absorption in the soft (0.2 to 2.4 keV) energy 
band, consistent with the mini-BAL nature of this quasar, while a comparison 
with ROSATPSPC observations from 1990 shows a drop of ~6.5 in the total soft 
bandpass flux. The observations also detect ~ 157 X-ray photons arising from 
extended emission of the nearby cluster (peaked ~ 42" SW of RX J0911. 4+0551) 
responsible for the large external shear present in the system. The Chandra ob- 
servation reveals the cluster emission to be complex and non-spherical, and yields 
a cluster temperature of kT = 2.3TJo'g keV and a 2.0 to 10 keV cluster luminosity 
within a 1 Mpc radius of L x = 7.6t° 6 2 x 10 43 ergs s 1 (error bars denote 90% 
confidence limits). Our mass estimate of the cluster within its virial radius is 
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2.3^ a y x 10 14 M©, and is a factor of two smaller than, although consistent with, 
previous mass estimates based on the observed cluster velocity dispersion. 

Subject headings: gravitational lensing: individual (RX J0911. 4+0551) — X-rays: 
general - quasars: individual (RX J0911. 4+0551) — quasars: (BAL) 

1. Introduction 

Gravitational lenses that produce multiple images of quasars can be used to measure 
cosmological parameters (Refsdal 1964; Kochanek 1996), to study the properties of the 
lensing galaxies (Keeton, Kochanek, & Falco 1998) and to resolve structures associated with 
the lensed quasar (Bunker, Moustakas, & Davis 2000). For the past twenty years such work 
has been carried out at optical and radio wavelengths. With the advent of the subarcsecond 
angular resolution of the Chandra X-ray Observatory (Weisskopf, O'Dell, & van SpeyBroeck 
1996), one can hope to carry out similar studies at X-ray wavelengths. As the X-ray emission 
from quasars has different spatial and temporal structure, X-ray studies offer a variety of 
new opportunities. In the present paper we report on Chandra observations of the quadruply 
imaged quasar RX J091 1.4+0551, which have implications both for the structure of the 
lensing potential and for the structure of the X-ray emitting region. 

To our knowledge, the observations of RX J0911. 4+0551 presented here are the first 
X-ray observations of the system since its detection by the RO SAT All Sky Survey in 1990 
October. The source was optically identified as a gravitationally lensed z = 2.8 quasar 
by Bade et al. (1997), and later confirmed as a quadruple system by Burud et al. (1998). 
Lensing models require a large external shear (7 > 0.15) to account for the system's complex 
image geometry. Burud et al. (1998) first suggested the source of this shear to be a group 
of galaxies in the southwest vicinity of RX J091 1.4+0551, which has been spectroscopically 
confirmed by Kneib et al. (2000) to be a high-redshift cluster of at least 24 galaxies at 
z = 0.769. The original R OSA T observation of RX J0911.4+0551 lacked the sensitivity and 
angular resolution required for a detailed X-ray analysis of the lensed quasar and cluster. 
Deconvolution of the Chandra observations presented here identifies all four components 
of the lensed quasar, and an adaptively smoothed image reveals the extent of the X-ray 
emission from the nearby cluster as well. The format of the paper is as follows: we describe 
the Chandra observations and reductions of RX J0911. 4+0551 in §2, discuss X-ray astrometry 
and relative photometry of the lens components in §2.1, and report spectral properties of 
the lensed quasar in §2.2. The X-ray detection and analysis of the cluster is presented in §3. 
We summarize our findings and briefly comment on possible sources of X-ray variability for 
the system in §4. 
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2. Observations and Analysis 

RX J0911. 4+0551 was observed for ~ 29 ks with the Chandra X-ray Observatory on 1999 
November 2. The data were obtained using the back-illuminated S3 chip of the Advanced 
CCD Imaging Spectrometer (ACIS; Garmire et al. 1992; Bautz et al. 1998) at a focal plane 
temperature of — 110°C. The telescope pointing placed RX J0911. 4+0551 at the default aim 
point for the ACIS-S array during the exposure. Three periods of background flaring were 
detected, during which time the count rate reached 2 to 4 times the nominal background rate. 
Since the contamination of point sources from background flaring is small, these times were 
not filtered from our observation except for the cluster analysis discussed in §3. Removal of 
the three X-ray flares reduces the effective exposure time of the data to ~ 26 ks. 

The data were reduced using the CIAO software tools, following the science threads 
outlined on the Chandra X-ray Observatory Center (CXC) user support homepage. The 
data used throughout this paper has been reprocessed by the CXC using software version 
R4CU5UPD9, which results in an improved aspect solution over earlier versions of the process- 
ing software. Also, in the standard CXC pipeline, the photon event positions are randomized 
by ± 0'.'246 (1 ACIS pixel = 0'.'492), which is mainly to reduce aliasing effects for observations 
of duration less than ~ 2 ks. Given the small angular extent of the lens (~ 3"), we repro- 
duced the event positions using the CXC tool acis_process_events without incorporating 
the randomization. 

Throughout this work, we have adopted a flat cosmological model parameterized by 
q Q = 0.5 and H = 50 km s -1 Mpc -1 . At the redshift of the lensing galaxy (z = 0.769; Kneib 
et al. 2000), this gives a physical scale of 8.2 kpc arcsec -1 . Unless otherwise noted, error 
bars denote 90% confidence limits. 



2.1. Relative Astrometry and Photometry 

In Figure la we show the X-ray image of RX J0911. 4+0551, resampled at a resolution of 
O'.'l pixel -1 and smoothed with a Gaussian of cr=0'.'2. A total of ~ 404 X-ray photons in the 
0.3 to 7.0 keV bandpass are estimated from the four lensed images of the quasar, which com- 
pensates for contributions from the X-ray background and from the nearby cluster. Adopting 
the nomenclature of Burud et al. (1998), component B is clearly resolved whereas compo- 
nents Al, A2 and A3 are not clearly separated. To enhance the image quality we applied 
the Richardson-Lucy maximum likelihood deconvolution technique (Richardson 1972; Lucy 
1974) to the X-ray image of RX J091 1.4+0551. The point spread function (PSF) incorpo- 
rated in the deconvolution was generated with the raytrace simulator MARX v3.0 (Wise 
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et al. 1997). The input spectrum for the simulation was based on the model of the ob- 
served spectrum of RX J091 1.4+0551 (see §2.2). In particular, we chose an absorbed power 
law model with Galactic absorption of N# = 0.036 x 10 22 cm -2 , intrinsic absorption at 
z = 2.8 of N# = 4.1 x 10 22 cm" 2 and a photon spectral index of 1.63. The PSF and image 
were then sampled at a resolution of O'.'l. In Figure lb we show the deconvolved image of 
RX J0911. 4+0551 smoothed with a Gaussian of a = 0'.'2. All four components are now re- 
solved with relative positions in good agreement (to better than / .'15) with those found from 
optical observations. Component B appears somewhat extended along the EW direction in 
both the raw and deconvoled images, which we attribute to the small number of counts (~ 
70 in total for 0.3 to 7.0 keV) detected for this component. In Table 1 we present the optical 
and X-ray offsets of the lensed images. 

The ACIS PSF (when properly normalized) gives the probability distribution of detect- 
ing a single X-ray photon at some position with respect to the center of a point source. One 
can therefore use several PSFs, of different relative intensities, to quantify the likelihood of 
obtaining the observed distribution of photon counts for a given set of image flux ratios and 
positions (e.g., see Cash 1979). Specifically, we have solved for the relative flux ratios of the 
system components using the Cash (1979) C statistic to quantify the relative goodness of fit 
between the observed and predicted number of photons per bin. The system was modeled 
using four of the MARX PSFs described above; relative positions of the four sources were 
held fixed according to the optical astrometry listed in Table 1, but the overall position of 
the system was allowed to vary during the fit. The best fit parameters were then determined 
by simultaneously varying the intensity of each component and the overall position of the 
system using Powell's direction-set method (Press et al. 1995) until a minimum C statistic 
was reached. This analysis was performed on the raw data at a resampling size of 0'.'06, or 
an eighth of an ACIS pixel. 

In Table 2 we present our optimal X-ray flux ratios for the 0.3 to 7.0 keV bandpass, along 
with contemporaneous (to within one day) /-band data obtained from an ongoing monitoring 
campaign of RX J091 1.4+0551 using the Nordic Optical Telescope (NOT; Hjorth et al. 2001). 
Since the separation between components Al and A2 is C'48, or slightly less than one ACIS 
pixel, it is prudent to discuss the combined flux of components Al and A2 when comparing 
to the optical data. Our maximum likelihood analysis yields an (Al + A2)/B X-ray flux 
ratio of 4.68 and an (Al + A2 + A3)/B flux ratio of 4.87, which are ~ 20% and ~ 30% 
smaller than the respective optical ratios of 5.67 and 6.95. A significant difference is observed 
for A3/B, where we find an X-ray flux ratio of 0.19 as compared to the optical ratio of 1.28. 
If we normalize the optical and X-ray fluxes at their respective (Al + A2) values, then 
component B is brighter by a factor of 1.2 in X-rays and component A3 is fainter by a factor 
of 5.7 in X-rays. We briefly discuss possible implications of the faintness of A3 in §4. 
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2.2. Spectral Properties 

A source spectrum of RX J0911. 4+0551 was extracted by combining events from all four 
quasar images and a background spectrum from events located within an annulus centered 
on RX J0911. 4+0551 with inner and outer radii of 15" and 30", respectively. Redistribution 
matrices and ancillary response files were generated using the CXC tools mkrmf and mkarf . A 
fit of an absorbed power-law model with intrinsic and Galactic absorption in the energy range 
0.5-6.0 keV yields a photon spectral index of I^IqJo an d a 0.2-2.4 keV flux (approximating 
the ROSAT bandpass) of fx = 3.31qi x 10~ 14 erg s _1 cm~ 2 . A fit of a simple absorbed power- 
law model with Galactic absorption of N# = 0.036 x 10 22 cm -2 , considering only energies 
above 1.5 keV, yields a photon spectral index of 2.22+Q4g, which is typical for high redshift 
radio-quiet quasars (Reimers et al. 1995). As we show in Figure 2, extrapolating this model 
to energies below 1.5 keV clearly indicates that considerable intrinsic absorption is present. 
This low-energy absorption is consistent with the mini-BAL nature of RX J091 1.4+0551 
(Bade et al. 1997), and has been observed in at least one other BAL quasar (e.g., PHL 5200, 
Mathur et al. 1995). For RX J0911.4+0551, the 0.2-2.4 keV flux of the unabsorbed model 
is fx = l-0loi x 10~ 13 erg s _1 cm~ 2 , and represents what we would expect to observe were 
the absorber not present. 

RX J091 1.4+0551 was observed with the ROSAT PSPC on 1990 October with a count- 
rate of 0.020 ± 0.0088 cnts s _1 (1 a). Assuming the simple absorbed power-law model 
from above, the corresponding X-ray flux for the ROSAT observation is fx = (2.2 ± 1.1) 
x 10~ 13 erg s~ x cm -2 (1 a). This includes emission from both the gravitational lens and 
extended emission from the nearby cluster (see §3). Within a circular aperture centered on 
RX J091 1.4+0551 of diameter V (twice the half-power diameter of the ROSAT PSPC), the 
Chandra estimate of the flux due to solely cluster emission is (6.6 ± 1.0) xl0~ 15 ergs cm -2 
s _1 (1 a). Therefore, the cluster contribution to the original ROSAT Qux was small, on the 
order of a few percent. When we compare the estimated ROSAT fiux from RX J0911. 4+0551 
to the Chandra observation, we find it has dropped by a factor of ~ 6.5 ± 3.3 (1 a) over the 
10-year period. Our spectral fits to the Chandra spectrum of RX J091 1.4+0551 indicate that 
a factor of 3.0 ± 1.5 (1 a) variation in the 0.2-2.4 keV flux is possible if one assumes that 
the presently observed intrinsic X-ray absorber were not present during the ROSAT PSPC 
observation. 



3. Cluster Analysis 

The Chandra observation of RX J091 1.4+0551 provides the first detailed X-ray study 
of the nearby cluster reported by Kneib et al. (2000). The extent of the X-ray emission from 



the cluster is apparent after adaptively smoothing the image with the CXC tool csmooth 
(Ebeling, White and Rangarajan, 2000). Contours for the smoothed image are displayed in 
Figure 3, and are shown overlaid on an UK composite exposure of the RX J091 1.4+0551 field 
from Burud et al. (1998). When creating the smoothed image, we have used a conservative 
lower bound on the signal to noise ratio (SN) of the smoothing kernel of SN^3, and have 
also excluded emission from within a small (~ 4" radius) region around RX J091 1.4+0551 
itself. The smoothed emission profile consists of a bright, extended component peaked ~ 42" 
southwest of component B, but is clearly elongated in the direction toward RX J0911. 4+0551. 
An iterative centroiding algorithm (Buote & Canizares 1992) applied to the binned event 
positions yields a centroid for the extended emission of Aq = — 12"0, AS = — 39"8, where 
offsets are with respect to the centroid of component B. The contour morphology suggests 
that the cluster mass distribution is complex and non-spherical, and possibly not dynamically 
relaxed. If the X-ray emission follows the mass distribution (which, simulations suggest, it 
need not always do during merging; see, e.g. Roettiger, Stone & Mushotzky 1998) then our 
image implies that the cluster is responsible for some fraction of the mass convergence at 
the location of the lens. 

Within an ellipse of 70" x 90" (dashed outline in Figure 3) encompassing the bulk of 
the cluster emission, we find a total of 157 net cluster counts above an estimated background 
of 222 counts. The corresponding 0.3 to 7.0 keV flux (assuming the spectral form described 
below) is f x = (2.5 ± 0.2) xlO" 14 erg s" 1 cm" 2 (1 a). 

To estimate an X-ray temperature for the cluster, we extracted a 0.3 to 7.0 keV source 
spectrum centered on the SW peak within a 20" aperture, and subtracted a background 
spectrum from an annulus with inner and outer radii of 90" and 140", respectively. Fitting 
the emission spectrum to an optically thin Raymond and Smith (1977) plasma model with a 
Galactic absorption of N H = 0.036 x 10 22 cm~ 2 , a metal abundance of 0.3 Z , and a cluster 
redshift of z = 0.769 (Kneib et al. 2000) yields a temperature of kT = 2.3+qJ at the 90% 
confidence level. This is consistent with the estimate of kT = 4.5 ±1.2 keV obtained by 
Kneib et al. (2000) from their measurement of the cluster velocity dispersion and the a — T x 
relationship (Girardi et al. 1996). 

We do not detect enough photons to constrain the radial surface brightness profile of 
the cluster. We obtain a crude estimate of the integrated cluster emission by assuming the 
profile is well-described by an isothermal (3— model (Cavaliere & Fusco-Femiano 1976) with 
parameters typical for nearby, non-cooling flow clusters {(3 = 0.6, r c = 0.22 Mpc; Mohr et 
al. 1999). Normalizing the model using the observed flux within the elliptical aperture, we 
predict a 2.0 to 10 keV luminosity within a 1 Mpc (122") radius of L x = 7.61™ x 10 43 ergs 
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From our temperature we may obtain a crude cluster mass estimate. We first estimate 
the virial radius for the cluster. Girardi et al. (1998) find, for nearby clusters, that the 
virial radius is Ry ~ 4<xiooo Mpc, where cr woo is the radial velocity dispersion divided 
by 1000 km s -1 . Assuming the virial radius scales with redshift as Ry ~ (1 + z)~ 3 / 2 , 
and given the observed velocity dispersion criooo — 0.84 (Kneib et al. 2000), we estimate 
Ry ~ 1.4/i5Q Mpc. Assuming an isothermal plasma with the density distribution described 
above and kT = 2.3 keV, the mass within the virial radius is M(R V ) = 2.3lJ;f x 10 14 M . 
Alternatively, the mass-temperature relation of Hojrth, Oukbir, & van Kampen (1998) gives 
M(R V ) = 3.2tH x 10 14 M . Here the 90% confidence limits are dominated by errors in the 
temperature, and do not include systematic errors arising from our ignorance of the slope of 
the surface brightness distribution, or from calibration uncertainties in the mass-temperature 
relationship. Our X-ray mass estimate is somewhat smaller than, although consistent with, 
the virial estimate of 4.8+2:0 x 10 14 M© (1 a; Q m = 1.0, Q A = 0.0) obtained by Kneib et 
al. (2000). Given the asymmetry of the X-ray emission, it is reasonable to speculate that 
the X-ray gas is not in hydrostatic equilibrium. In this case, one might well expect the 
emission weighted X-ray temperature to lead to an underestimate of the gravitating mass 
(e.g. Roettiger, Stone and Mushotzky 1998). In any event, given the large uncertainties, the 
agreement among the mass estimates must be regarded as good. 



4. Discussion and Conclusions 

Our analysis has detected two indications of X-ray variability for RX J0911. 4+0551. 
First, component A3 is dimmer in X-rays by a factor of ~6 when compared to contem- 
poraneous optical data. The reason for A3's X-ray faintness is not known. One possible 
explanation is differential N H column densities along the lines of sight to the quasar images. 
Although an attempt was made to search for differential X-ray extinction between the four 
lens components (the details of which will be presented in a forthcoming paper; Chartas et 
al. 2001), a preliminary analysis is inconclusive; the intrinsic Nh column density of image 
A3 is poorly constrained by the Chandra data. A demagnifying microlensing event of com- 
ponent A3 is also a possible, although somewhat unlikely, explanation. A six-fold drop in 
flux would require fluctuations of the order Am ~ 2 magnitudes; microlensing dips of this 
size are rare in theory (Witt, Mao, & Schechter 1995). 

Second, the total observed flux for the system has dropped by a factor ~ 6.5 when 
compared to ROSAT PSPC observations obtained in 1990 October. Such a drop in flux 
may be a result of quasar engine variability. Almaini et al. (2000) have observed engine 
variability of a factor of three in flux over a period of several days, so a factor of ~ 6.5 
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between the RO 'SAT and Chandra observations is not unlikely. 

In conclusion, the unique imaging capabilities of Chandra has resolved all four images of 
the quadruply lensed system RX J0911.4+0551 (detecting ~404 photons in the 0.3 to 7.0 keV 
energy range), with relative image positions in good agreement with optical measurements. 
In addition to the two indications of X-ray variability described above, we have also detected 
strong soft (< 1.5 keV) X-ray absorption in the spectra of the combined quasar images, 
which is consistent with the mini-BAL nature of the quasar. The Chandra observation 
of RX J0911. 4+0551 has also identified extended X-ray emission from the nearby cluster, 
detecting ~157 photons (0.3 to 7.0 keV) in a region of ~ 1.4 arcmin 2 . The morphology of 
the cluster emission suggests a complex and non-spherical cluster mass distribution. Our 
estimates of the cluster X-ray temperature (kT = 2.31q| keV) and virial mass (2.3^7 x 
10 14 M ) are consistent with previous estimates derived from the observed velocity dispersion 
of the cluster galaxies (Kneib et al. 2000). 
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Fig. 1. — Panel a: Chandra X-ray image of the quadruply lensed quasar RX J0911. 4+0551, 
resampled at one-fifth of an ACIS pixel (1 ACIS pixel = 0'.'492) and smoothed with a Gaussian 
of a=0'!2. Panel b: Deconvolution of the Chandra X-ray image using the Richardson-Lucy 
algorithm. The deconvolved image has been smoothed with a Gaussian of a=0'/2. The four 
components of the system are labeled according to the nomenclature of Burud et al. (1998). 
North is up and East is left; the width of each panel is ~ 5.5". See inset of Figure 3 to 
compare the image configuration with an HST/NICMOS observation of the system. 

Fig. 2. — Top panel: Chandra spectral data of the combined quasar images of 
RX J0911. 4+0551 with a best-fit power law (including Galactic and intrinsic absorption) 
to the hard (1.5-7.0 keV) energy band. Bottom panel: Deviation (in units of one standard 
deviation) of the observed spectral counts to the hard energy band fit. 

Fig. 3. — Solid lines: X-ray contours (logarithmically spaced) of an adaptively smoothed 
Chandra image of the RX J0911. 4+0551 field (excluding emission from RX J0911. 4+0551 
itself). The background image is a composite optical image taken with the NOT and 
ESO/NTT (Burud et al. 1998). X-ray contours have been aligned with the optical ex- 
posure according to the Chandra aspect solution. Dashed line: Elliptical aperture used in §3 
to estimate the cluster luminosity. Inset: Closeup of RX J091 1.4+0551 taken from archival 
HST/NICMOS imaging of E. Falco (PI). The 4 quasar images and lensing galaxy are labeled. 



Table 1. Optical and X-ray Offsets of RX J0911.4+0551 Components 



Telescope 


B 


Al 


A2 


A3 




Aa("), A5 (") 


Aa("), A5 (") 


Aa("), A5 (") 


Aa("), A5 (") 


NOT 


0.00, 0.00 


2.935±0.004, -0.785±0.009 


3.194±0.007, -0.383±0.007 


2.922±0.008, 0.161±0.009 


Chandra 


0.00, 0.00 


2.88 ±0.20 , -0.88 ±0.20 


3.32 ±0.20 , -0.41 ±0.20 


2.88 ±0.20 , 0.17 ±0.20 



Note. - Optical positions and corresponding 1 a error bars are from observations obtained with the Nordic 
Optical Telescope (NOT; Burud et al. 1998). Relative offsets are in the sense of component A minus component B. 
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Table 2. Contemporaneous Optical and X-ray Flux Ratios of RX J0911. 4+0551 

Components 



Waveband 


(Al+A2)/B 


(Al+A2+A3)/B 


(A3)/B 


I band 
0.3-7.0 keV 


5.67 ± 0.14 
4 68 +0 - 30 


6.95 ± 0.16 

4.87 ±g : g 


1.28 ± 0.07 

0.19 ±8# 



Note. - - /-band ratios and corresponding la error bars are 
from an ongoing monitoring campaign of RX J0911. 4+0551 
with the Nordic Optical Telescope (Hjorth et al. 2001). Error 
bars for the X-ray ratios represent 90% confidence intervals. 
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